Hard-core Bose-Fermi mixture in one-dimensional split traps 
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We consider a strongly interacting one-dimensional (ID) Bose-Fermi mixture confined in a hard 
wall trap or a harmonic oscillator trap with a tunable 5-function barrier at the trap center. The 
mixture consists of ID Bose gas with repulsive interactions and of ID noninteracting spin-aligned 
Fermi gas, both species interacting through hard-core interactions. Using a generalized Bose-Fermi 
mapping, we calculated the reduced single-particle density matrix and the momentum distribution 
of the gas as a function of barrier strength and the parity of particle number. The secondary peaks 
in the momentum distribution show remarkable correlation between particles on the two sides of 
the split. 
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I. INTRODUCTION 

In recent years the strongly interacting ultra-cold 
atoms system in quasi-one dimensional(lD) have at- 
tracted great interests both in experiments and theories. 
With two perpendicular optical lattices to realize the 
ID system and Feshbach resonance or the confinement- 
induced resonance to tune the ID effective interact- 
ing strength, Tonks-Girarideau(TG) gas [H, [1], a Bose- 
Einstein condensate in which the repulsive interactions 
between bosonic particles dominate the physics of the 
system, has been realized in experiments. For strong at- 
tractive interactions between atoms, the Super Tonks- 
Girardeau gas [3] represents an excited quantum gas 
phase in a ID spatial geometry. Sudden switching from 
infinitely strong repulsive to infinitely attractive interac- 
tions stabilizes the gas against collapse and connects the 
ground state of the Tonks gas to the excited state of the 
Super Tonks gas. Meanwhile, the study of a two spin 
mixture of ID strongly interacting fermions gas at finite 
spin imbalance demonstrates how ultracold atomic gases 
in ID may be used to create non-trivial new phases of 
matter experimentally [1, Q . Despite intense theoretical 
and experimental efforts, however, the strongly interact- 
ing Bose-Fermi mixtures remains largely elusive, where 
interesting phenomena such as phase separation are pre- 
dicted to occur [^-Q. 

Most of the theoretical research on Bose-Fermi mix- 
tures so far has been concentrated on three dimensional 
systems, and only recently ID systems started attract- 
ing attention. Theoretical investigations on the quasi-lD 
Bose-Fermi mixture so far have focused on the phase dia- 
grams, ground-state and thermodynamical properties in 
the scheme of Luttinger liquid theory [HI and Bethe 
ansatz method Il2j |. For the ID ultra-cold atoms sys- 
tem with infinity repulsive interaction between particles, 
the Fermi-Bose Mapping (FBM) theorem is used to de- 



rive the wavefunctions and ground state properties of 
the system. In the simplest case the mapping function 
relates the system of impenetrable bosons with that of 
non-interacting fermions jl3| and similar schemes have 
been extended to various hard-core systems like spinor 
Bose gases spin- 1/2 fermions [3| and Bose-Fermi 
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mixture 

The starting point of the many-particle solutions of a 
given geometry for the TG gas is the exact single-particle 
eigenstates. Given the limited examples of exactly solv- 
able single-particle problems in quantum mechanics, ex- 
actly solvable many-particle problems are also rare even 
in the TG limit. Here we choose two typical split traps 
to load the mixture, i.e., a split hard wall and a split 
harmonic oscillator. Investigating the model of the 5- 
split potential can be justified in several ways and related 
works have been done jl7l - [20j . For mesoscopic systems, 
the ground state properties qualitatively differ for sys- 
tem with a split at the trap center, especially when the 
number of particles is odd. On the other hand, (5-split 
potential may be viewed as a generic model for double- 
well structures or, alternatively, as a good approxima- 
tion to the problem of a trap with an impurity at the 
center. Experimentally this kind of split is easy to real- 
ize by adding an additional laser in the trap center pH . 
Tunneling of particles through an energetically forbidden 
region reveals more and more experimentally accessible 
examples of nontrivial quantum phases [22-25]. 

The paper is organized as follows. After a brief intro- 
duction of the many-body system, in Sec. II we review 
the Fermi-Bose mapping theorem for Bose-Fermi mixture 
and describe the associated single-particle eigenfunctions 
and eigenvalues of the ^-split traps. In Sec. Ill we cal- 
culate the reduced single-particle density matrix and in- 
vestigate the influence of the splitting strength on the 
momentum distribution which is experimentally acces- 
sible. We observe the emergence of bimodal secondary 
peaks at the neck of the central peaks for various barrier 
strengths and analyze the parity effect of the total num- 
ber of atoms. Finally, in Sec. IV we make concluding 
remarks. 



2 



II. BOSONS AND FERMIONS IN SPLIT TRAPS 

Consider a mixture of Nb bosons and Np fermions 
trapped in a tight atomic waveguide which restricts the 
dynamics of the system into a quasi one-dimensional sys- 
tem in the longitudinal direction. Assume that the in- 
teraction potentials are short-ranged, the many-particle 
Hamiltonian at low linear density can be written as 

+ ^3bb*b*I*6*b + 5fc/^i*/*/*b|> (1) 

Here, 5*6, 5*/ are the boson and fermion field operators. 
The bosonic and fermionic particles are assumed to share 
the same mass mi, = nif = m, which could be realized 
by choosing different isotopes of a given alkali element, 

e.g. 40^_39(41)^ ^ 86(84)^^ _ 87(85)^^ ^ 

model, the particles experience both boson-boson (gtb) 
and boson-fermion (gbf) contact interactions, while the 
fermion-fermion interaction is forbidden by the Pauli ex- 
clusive principle {gjf = 0). We also assume the trap- 
ping potential are exactly the same for both fermions 
and bosons, Vb{x) = Vf{x) = V{x). 

A. The Fermi-Bose mapping 

The original Fermi-Bose mapping theorem only related 
strongly interacting bosons to ideal fermions [l3| . It 
was recently found that the mapping idea can be ap- 
plied to Fermi Tonks gas and to a series mixture sys- 
tem as well [l^. Concentrating on the situation rele- 
vant to our system, we denote the space coordinates by 
X = (xi, X2, xjv), where Xi is the coordinate of a boson 
when i g {1, Nb} , or else it labels a fermion provided 
that i e {Nb + 1, Nb + Np = N}. Under the hard- 
core condition gbb,gbf — >■ oo, we may safely drop the two 
interaction energies in the Hamiltonian and as a result 
the many-body wave function of the system is subject to 
a constraint that 

$(X) = 0, if X, = Xj (2) 

for i ^ j. The first quantized Hamiltonian of our bose- 
fermi mixture is then a sum of one-body operators, H — 
Yl!i=ihi, where 

h = -^dl+V[x,). (3) 

Using the Fermi-Bose mapping theorem [l^, the to- 
tal wave function satisfying the constraint is con- 
structed as 

MX)=A{X)<i>D{X), (4) 



where <^d{X) is a Slater determinant of N orthonormal 
orbitals U2{x), . . . WAr(x) occupied by N particles 

<^d{X) = Y,^{P)ui[P^i)MPx2) . ■ -UNiPxN). (5) 
p 

The sum runs over all iV! possible permutations of these 
variables including permutations exchanging bosons with 
fermions, and s{P) — ±1 for even/odd times of permu- 
tation. The symmetry of the system is repaired by the 
mapping function A{X) 

Nb N 

MX) = W sgn(xj-xi)'^ ]^ sgn{xj - xi) 

1<]<1<Nb 3 = 1 1=Nb + 1 

(6) 

so that in general the wavefunction is symmetric under 
permutations of bosons and antisymmetric under permu- 
tation of fermions. The sign function sgn[x) is -|-1(— 1) 
for X > 0(x < 0). 

B. Eigenstates and eigenvalues of the (J-spHt traps 

The orthonormal orbitals Un{x) are eigenfuctions of 
the single-particle Hamiltonian with eigenvalues tn 

(^-^^x + V (x)^ Ur,{x) = e„M„(x) (7) 

with n = 1,2, ... N . In this paper we consider the case 
that the mixture is subjected to two kinds of split ex- 
ternal trapping potentials, that is, (A) a hard wall trap 
which is zero in the region (—a, a) and infinite outside, 
and (B) a harmonic oscillator potential with frequency 
Lo, both with a J-type barrier located at the origin a; = 0. 
We denote the strength of this barrier by a positive pa- 
rameter K. Thus 

for model (A) and 

V (x) = \nij?x^ -f K(5 (x) (9) 

for model (B). The single-particle eigenstates of the delta- 
split hard wall trap can be found in quantum mechanics 
textbook, and those of the delta-split harmonic oscillator 
have recently been discussed in ref 28j and we will briefly 
review the solutions here for completeness. 

The eigenfunctions are either symmetric or antisym- 
metric due to the parity of both potentials. For model 
(A), the analytic symmetric eigenfunctions are 

Un {x) = C [ cos {kx) -\ — — s\n{k \x\) I (10) 

for |a:| < a, and it„ {x) = for \x\ > a. Here, C is the nor- 
malization constant and k is the wave vector of the par- 
ticle, determined by tan (/ca) = —h^k/rriK. The eigenen- 
ergies E = fi.^fc^/2m are given by the graphical solutions 
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of k, that is, the intersection points of the straight line 
—h^k/mn and the tangent function which can be labeled 
as n= 1,3,5.... Ask— >-0, the eigenenergies reduce 
to those of the ordinary infinite well of width 2a. As 
K — > oo, the barrier becomes impenetrable, and we have 
two isolated infinite square wells of width a. By contrast, 
the antisymmetric eigenfunctions are zero at the origin, 
so they never "feel" the delta function at all. The eigen- 
functions and eigenenergies are exactly the same as those 
for the infinite square well of width 2a 



Un (x) = C sin {kx) 



(11) 



with k ~ n-njla and En — n^Ti'^h^ /2m{2aY {n — 
2, 4, 6, . . .). We will use dimensional variables in units of 
the half split square potential length a, all momenta in 
units of p = ?i/2a, and all energies in units of 7r^?i^/2ma^. 

For model (B), we immediately know that the antisym- 
metric eigenfunctions of the simple harmonic oscillator 
remain good eigenfunctions for the 5-split oscillator, as 
they vanish at the exact position of the disturbance. We 
have 



where Nn 

Q — X I Xq 



{^2^n\xoscy^'^ 



(12) 



with n = 1, 3, 5, . . . Here 



= ^/h/muj and Hn{Q) are the n-th 
order Hermite polynomials. The corresponding energies 
are given by the eigenvalues of the odd parity states of 
the harmonic oscillator En = {n + 1/2) Tllu. 

The even eigenstates of the simple harmonic oscillator, 
on the other hand, have an extremum at a; = 0. They are 
therefore strongly influenced by the splitting potential 
and can be found to be [H, [l!] 

(13) 

where U{a, b, z) are the Kummer's confluent hypergeo- 
metric functions [l^ and Vn = En/2huj — 1/4 is the 
non-integer analog of the principle quantum number n 
of the harmonic oscillator. The eigenenergies En {n = 
0, 2, 4, . . .) are given by the solution of 



Unix) - Ar„e-Q'/2f^(-M„, 1/2, Q2) 



2T{~Vn + 1/2) + KT{~yn) = 



(14) 



Similarly in the following the variables are rescaled 
in units of the harmonic-oscillator length Xosa the 
harmonic-oscillator momentum posc = Ti/xosc = Vmhtu, 
the harmonic-oscillator energy hio. 

With these single-particle eigenstates, we are now in 
a position to build the Slater determinant $£)(X) for 
a system of noninteracting fermions. The Fermi-Bose 
mapping theorem then allows us for the calculation of 
the exact many particle wave function ^(X) from the 
fermionic result. 



III. MOMENTUM DISTRIBUTIONS OF THE 
BOSE-FERMI MIXTURE 

The quantum correlations of the ID quantum gases 
can be obtained from the reduced single-particle density 
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FIG. 1: (Color online) Total density distributions for odd 
number of atoms (A*' = Nb + Nf ~ 5) at different strengthes 
of the split barrier (a) for Model A (x in units of a), (b) for 
Model B [x in units of ^jTilmbS). 



(RSPD) matrices defined as 

Pb{x.V) = NbJ dX'<i>*{x,X')<i>*{y,X') 
PFix,y) ^ Nf I dX"^*(X",xmX'\y) (15) 



for bosons and fermions, respectively. Here X' = 
(x2, xat) and X" = (xi, ...,Xn-i)- In the TG inter- 
action limit, however, the diagonal elements are nothing 
but the single-particle density profiles pb{x,x) — Pb{x) 
or pp{x,x) = Pf{x) which are exactly the same up to 
normalization factors due to the fact that the square of 
the mapping function A{X) in (|6]) is unity. They both are 
proportional to the spatial density ptg{x) of a trapped 
TG gas made of N bosons. In Fig 1 we illustrate the 
normalized density profiles 



Pb{x)/Nb = Pf{x)/Nf = Ptg{x)/N 



(16) 



for split hard wall trap (Model A) and split harmonic os- 
cillator (Model B), respectively. Here we are interested 
in the effect of the ^-barrier in the center and results 
for four values of the barrier strength k, ~ 0.1, 1, 10, 100 
are presented. It has been shown that for even number 
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of particles the barrier separates the particles into two 
half wells and increasing the barrier strength leads to the 
emergence of a quadrant separation of which the inter- 
ference is negligible [1, d^. Thus only results for odd 
number of particles are given here and we found that 
the density profile relies on the total number of parti- 
cles, instead of the number of bosons and fermions. For 
negligibly small split barrier, the density profiles present 
exactly the same number of density maxima as the total 
number of atoms. The corresponding density profile in 
the large- iV limit is flat in the square well and is a semicir- 
cle of radius R = y/2N in the harmonic trap. Increasing 
the barrier strength would gradually cut the density pro- 
files into two symmetric parts and the number of peaks 
do not match the number of particles any more. This 
is a clear signature of the coherence inherent in the sys- 
tem, which becomes more pronounced in the momentum 
distribution. The delta-splitting seems more efficient in 
suppressing the density in a harmonic oscillator. This can 
be understood by noticing that the energy unit adopted 
in hard wall is 7r^/2 times larger than that in harmonic 
oscillator. 

The Fermi-Bose mapping theorem gives identical den- 
sity distributions of a sample for bosons and fermions in 
the TG limit. However, the momentum distribution can 
still be used for distinction, which can be calculated from 
the reduced single-particle density matrix 



1 

2^ 



Pb,f 



{x,y)e-'P^''-y^dxdy (17) 



which is normalized to N. The momentum distribution 
of the bosonic component nsip) is equivalent to that 
of N TG bosons placed within the split traps. This 
again comes from the fact that the symmetry of both 
boson-boson and boson-fermion permutations has been 
repaired in the mapping function and has already been 
addressed in Ref. [8l lla, Hi]. The physical reason for 
this identical momentum distribution originates from the 
equal- weighted superposition of the orthonormal orbitals 
in the construction of the many-body wavefunction ([5]). 
Thus we now turn our attention to the behavior of the 
fermionic momentum distribution in the mixture np{p) 
as a function of the split barrier strengths. The deter- 
mination of the fermionic momentum distribution proves 
very complex even in the homogeneous mixture and the 
extension to the harmonically trapped mixture has been 
done in Ref. [S^l. The calculation is more involved be- 
cause the expression for the density matrix could not 
be reduced to a simpler analytical formula, and we have 
resorted to numerical computations. Following a more 
convenient scheme, we use the alternative expression for 
momentum distributuion 



npip) ^ N dX" ^{X",p) 



(18) 



Here ^{X",p) is the Fourier transformation of 




FIG. 2: (Color online) Momentum distributions for even num- 
ber of atoms (A'^ = Nb + Nf = 4) at strong split barrier 
K = 100 (a) for Model A (p in units of h/2a), (b) for Model 
B [p in units of \/ mhco). 



^{X" , xn) with respect to the last fermionic variable xn 

^ /• + 00 



/27r 



dxN<^{X",XN)e-'P''". (19) 



In a recent paper, Lelas et al. studied the RSPDM, mo- 
mentum distribution, natural orbitals and their occupan- 
cies for a Bose- Fermi mixture in an alternative form of a 
harmonic potential with a Gaussian-type barrier in the 
center, however, for a rather strong fixed barrier, i.e. the 
barrier is at least several times larger than the energy 
of the iV-th single-particle state of the potential. We 
present here the results for various barrier strengths and 
analyze the dependence of the distribution profiles on 
the parity of the total number of particles. Furthermore 
our results provide a direct comparison of the momen- 
tum distribution for different type of trapping potentials 
(hard- wall and harmonic oscillator). 

Figure 2 shows our results for a mixture with a fixed 
even total number of iV = 4 particles when the number 
of bosons is increased from Nb = up to 3. A mixture 
with Nb = N 01 N — 1 makes no difference because one 
can not imprint any statistical property onto the only left 
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FIG. 3: (Color online) Momentum distributions for odd num- 
ber of atoms (A'' = Nb + Nf = 5) at weak split barrier k = 0.1 
(a) for Model A (p in units of h/2a), (b) for Model B (p in 
units of Vmhuj). 



FIG. 4: (Color online) Momentum distributions for odd num- 
ber of atoms (A^ = Nb + Np = 5) at strong split barrier 
K = 100 (a) for Model A (p in units of h/2a), (b) for Model 
B (p in units of V mhoj). 



atom. Both distributions for split hard wall and split har- 
monic oscillator take smooth bell-shaped profiles, with 
explicit difference being the characteristic momentum up 
to which the fermionic momentum distribution is signif- 
icantly different from zero. The signature of fermioniza- 
tion becomes more evident for less bosons. We also note 
that the split barrier in the center hardly affect the mo- 
mentum distribution of even number of atoms, except a 
little broadening of the momentum distribution and low- 
ering of its peak. 

Things become much more different for odd number of 
particles. We first illustrate the situation for a weak bar- 
rier. Figure 3 shows our results for a BF mixture with a 
fixed total number of A'^ = 5 particles when the number 
of bosons is increased from Nb = up to 4 (or, equiva- 
lently, 5). We observe Np "fermionic" oscillations devel- 
oped in the momentum distribution, yet the oscillations 
in harmonic oscillator are more prominent. Fermionic os- 
cillations are suppressed for more bosons and the tails of 
the distribution become more easily seen. Finally, when 
the atoms in the mixture are all bosons, we see typical 
Tonks-Girardeau-type momentum distribution. 



A strong split barrier greatly modifies the momentum 
distribution of odd number of atoms. Quite similar to 
the bosonic case, one can see the emergence of bimodal 
secondary peaks at the neck of the central peaks in Figure 
4, which stems from the interference of the particles in 
two almost separate wells. The secondary peaks remain 
prominent for all combination of bosons and fermions, 
moving inward when the number of bosons is increased. 
Again we observe that the characteristic momentum is 
notably different for the square well and for the harmonic 
confinement. For the square well it is of the order of fc ^ 
ttN/2, while for the harmonic trap it goes as fc ^ yj2N . 
The large- asymptotic in absence of bosons in both 
cases is step function. 



IV. CONCLUSION 

In conclusion, we have shown the ground state prop- 
erties of a Bose-Fermi mixture in split potential wells. 
We observe that the most significant difference between 
the total even and odd numbers of particles occurs in the 
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emergence of bimodal secondary peaks in the momen- 
tum distributions. Three critical features can be seen 
(1) the momentum distribution depends on the ratio of 
the bosons and fermions in the mixture, as well as the 
total number of the atoms; (2) the presence of a strong 
split barrier at the trap center enhance greatly the cor- 
relation of the atoms on the two separated wells; (3) the 
dependence of characteristic momentum on the number 
of atoms relies on the shape of trapping potential. 
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